Технические науки 


УДК 541.138.2 

МОДЕЛИРОВАНИЕ ЭЛЕКТРОДНЫХ ПРОЦЕССОВ ПРИ ПРОТЕКАНИИ ПРОСТОЙ 
ЭЛЕКТРОХИМИЧЕСКОЙ РЕАКЦИИ И ЛИНЕЙНОМ ИЗМЕНЕНИИ ПОТЕНЦИАЛА 
НА ЭЛЕКТРОДАХ ОГРАНИЧЕННОГО ОБЪЁМА 

И. С. Антипенко, В. И. Косинцев, Ю.А. Иванов*, Б.Ф. Назаров* 

Томский политехнический университет 
*НПП «Техноаналит», г. Томск 
Е-таіІ: ІаѲтаіІТотзкпеІги 

Показана возможность расчёта вольтамперных кривых при обратимом электродном процессе на плоских электродах для любых 
значений безразмерных параметров Н и Ѳ 0 , характеризующих природу электродного процесса, параметры электрода, скорость 
и форму изменения потенциала за время, соизмеримое с реальным временем проведения эксперимента. Величина Н показы- 
вает влияние соотношения скорости электрохимической реакции, определяемой скоростью изменения потенциала, диффузии, 
обусловленной взаимодействием между молекулами в растворах (металлических, жидких) и квадратом толщины плёнки элек- 
трода (анодный процесс) или электролита (катодный процесс). Величина Ѳ 0 показывает влияние соотношения скоростей элек- 
трохимических реакций в зависимости от равновесного потенциала. 


Имеется много работ, посвященных изучению 
обратимых электродных процессов при линейном 
меняющемся потенциале на плоских электродах 
ограниченного объёма. Наибольший вклад в изуче- 
ние этого вопроса внесли авторы работ [1-6], в ко- 
торых проведены расчеты различных уравнений для 
системы уравнений в частных производных в усло- 
виях ограниченной-полубесконечной диффузии. К 
недостаткам этих работ следует отнести громозд- 
кость вычислений, которые невозможно проводить 
на современных персональных компьютерах. 

Целью данной работы является моделирование 
обратимого электродного процесса при протека- 
нии простых электродных реакций и линейном из- 
менении потенциала для любых значений параме- 

С° пР 

тров Н=(пР1 2 1Ѵ)/(КТВ) и ѳ 0 = = ехр (Е р -Е°) — 

С к К1 

- отношение начальных концентраций окисленной 
(СЗх) и восстановленной (С а 0 ) форм вещества, а так- 
же максимальное упрощение расчётов полных воль- 
тамперных кривых при заданном значении равно- 
весного потенциала. Все величины, входя щ ие в па- 
раметр Н и 0 О , являются разнородными (/-толщина 
плёнки электрода, п - число электронов, прини- 
мающих участие в электродном процессе, ]Ѵ- ско- 
рость изменения потенциала, О - коэффициент 
диффузии, Е - стандартный потенциал, Е р - равно- 
весный потенциал, Т- температура, Р- постоянная 
Фарадея, К — универсальная газовая постоянная), 
поэтому с точки зрения процессов и аппаратов хи- 
мической технологии параметры Н и 0 О являются 
критерием подобия электрохимического процесса. 

Для линейной ограниченной-полубесконечной 
диффузии справедлива следующая математическая 
модель [1—4]: 

дС к (х,і) _ д 2 С к (х,і ) 

ді к дх 2 

с начальным 1=0: 

С к ( х , 0) = = соші 


( 1 ) 

( 2 ) 


и граничными условиями /> 0: 

дС к (х,і)\ 


С(1,1) = С'І = ./(/), (4) 

где х=0 - положение внутренней границы элек- 
трод-подложка; х=1 - положение границы элек- 
трод - раствор. Начальное условие (2) соответству- 
ет равномерному распределению концентрации 
металла в плёнке перед началом его электрораство- 
рения. Граничное условие (3) указывает на отсут- 
ствие потока вещества через границу раздела фаз 
(плёнка - подложка электрода). Условие (4) отра- 
жает произвольность концентрации вещества 
(С к (/,/)=Ся) на поверхности электрода, что эквива- 
лентно произвольности задания потенциала (т. е. 
значение потенциала меньше того, при котором 
достигается предельный диффузионный ток). Ре- 
шая ур. (1) с условиями (2-4) методом операцион- 
ного исчисления (преобразования Лапласа), полу- 
чим в пространстве изображений следующие выра- 
жения для концентрации вещества в любой точке 
электрода (ртутной пленки): 


С к (х, $) =—+ ($ • /($ ) - С” ) 

5 



и для потока вещества через границу электрод - 
раствор: 




<іС к (х,$) 


дх 


I х=1 



где У - площадь поверхности электрода, 5 - вспо- 
могательная переменная. 
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Переходя из пространства изображений в про- 
странство оригиналов, получим следующее выра- 
жение для потока: 


Б к 8 


с!С к (/, і) 
сіх 



дт 


9(1,2)1іт 
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( 6 ) 


где — кѵ - оригинал от выражения, стоящего в 

дт 

скобках в уравнении (5), а д Мѵ ((—т) - оригинал от 
остальной части выражения (5), причем он соот- 
ветствует предельному потоку на пленочный элек- 
трод с учетом замены ?на (/—г). Для того, чтобы пе- 
рейти от потока к току, умножим выражение (6) на 
величину пР, после чего получим: 

пМуІЩ-С^ = 

дх ^ 

= пРБ с 1\^- Чмш {1-т)ат, (7) 


і = пР8 
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дх 
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где С'и~ - нормированная концентрация ме- 

Оі 

талла на электроде. 

Значение определяется формулами: 
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I к= 0 ^ * 


г° 

г 8 - 

'-'О 


і+ѳ 


і +Ѳо -М^) х 

0 і+ѳ 

х ХІ — 1* ' 

&{і + Ѳ 


2 і/ Л 

к-Н /г 


{аіу 


2 2 


( 11 ) 


Из выражения (11) можно вычислить производ- 
ную: 

дС^(1,аІ) _ С К (/, сг2) — С к (/, а(1 + А()) 

д(аі) д '->° А (а() 


подставив которую в формулу (7), можно рассчи- 
тать зависимость тока от времени (потенциала). 
Аналитическое выражение для производной не 
приводится из-за слишком громоздкого её вида. 

Выражение для анодного тока в компактном 
виде будет иметь вид: 


і = 


п 2 Р 2 8С° к ПѴ . . 
к Х(а*), 
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- безразмерная функция тока. 


( 12 ) 


(13) 


или ^^ 2)Ш 



а + 2Х(-1)*ехр 

к = 1 


2 


V 


к 2 ГА 




(9) 


Заметим, что формулы (8) и (9) тождественны, 
причем формула (9) удобна для расчетов при боль- 
ших, а (8) - при малых значениях / (или Н). 

В выражении (5) величину д*"а( г ) необходимо 

дт 

доопределить, исходя из дополнительных краевых 
условий и уравнения Нернста. 


В работах [2,5] получено дополнительное крае- 
вое условие для ограниченной-полубесконечной 
диффузии, которое имеет вид: 


Рр 2 С'о + Рі 2 С К 

оАс° К 


1 + Ѳ 0 
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егіс 
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где: ѳ = % = ехр( (Е р - Е° + Ш) —1 - модифици- 

руя ѵ РТ ) 

рованное уравнение Нернста. 

пЕ 

Обозначим суі = 1 = — (Е-Е ) как без- 

КТ КТ р 

размерный потенциал. 

Решая ур. (10) совместно с уравнением Нерн- 
ста, получим 


Уравнение тока для катодного процесса с уче- 
том соотношения Сд=Ѳ 0 СІ и знака изменения по- 
тенциала имеет вид х((ті)=х((7І)/Ѳ І] . 

Для плёночного электрода справедливо: 


Яш (И’ 0 


2 
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где Р„=(2«-1)— , то очевидно, что при Я->0 
Нт с /ііт (2Я,/) = о, при ЯО и \іт ц ш {Н,і) = <х> при 

#-»0 Я-> о 

1= 0, а также 



-и;,— і 

" Я 


т = 1 . 


То есть, д м (Н,1) при Н -> 0 является дельтаобраз- 
ной последовательностью по определению [7]. С 
учетом свойств дельта-функции [8] выражение (13) 
примет вид: 


Пт у(оД) 


сі(сті) 


(14) 


а выражение (12): 
п 2 Р 2 8СІПУ 


і = ■ 


КТ 




п 2 Р 2 8СІт ( ІС(1,оі ) 
КТ с/(аІ) 


( 15 ) 


Отметим, что формула (13) впервые была полу- 
чена В. С. Баевым [9] без детального анализа. 
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Естественно, что 8СЦ1=(), поэтому соотноше- рис. 1-3 представлены результаты расчётов безраз- 
ние (15) можно записать в виде: мерной функции тока ^(сті) в зависимости от вели- 

п 2 Р 2 уу сІС'(І а() чины безразмерного потенциала с ті, рассчитанные 

і = — — — ( 2 — - — — — • (16) для различных значений Н и 0 О . 

Исследованы форма ВАК в зависимости от Н (от 
Дальнейшие расчёты проводились для значения о, 0004 до 1 ) для анодного и катодного процесса с до- 
0= 1 , то есть при постоянной площади под вольтам- полнительным учетом влияния равновесного потен- 
перной кривой (ВАК), которая не зависит от пара- циала . Проведённые расчёты по формулам (14, 15) и 

метров Н и Ѳ п . Таким образом, нами исследовано (12, 13) показывают, что уже в реальных условиях ме- 

влияние граничных условий (параметр Н) и началъ- хода инверсионной вольтамперометрии пользоваться 
ных условий (0 О ) на форму и подобие ВАК. На формулами (14, 15) предпочтительнее, чем точными 



О 2 4 6 8 10 

Безразмерный потенциал. оГ 

Рис. 1. Расчётные вольтамперные кривые по формуле ( 14) для различных значений параметра Н 



О 2 4 6 8 10 

Безразмерный потенциал, оі 


Рис. 2. Относительная разница между вольтамперными кривыми, рассчитанных по формуле (14) и ( 13) при варьировании па- 
раметра Н 



0 2 4 6 


Безразмерный потенциал, оі 

Рис. 3. Влияние начального равновесного потенциала на форму вольтамперных кривых, рассчитанных по формуле (16) для 
различных значений параметра Ѳ 0 
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формулами (12) и (13). Максимальная ошибка вычи- 
слений по формулам (14), (15) не превышает 15 % и 
связана со смещением расчётной ВАК на 3... 5 мВ в от- 
рицательную сторону. Скорость расчета возрастает на 
два-три порядка. При совмещении вольтамперных 
кривых по оси потенциалов ошибки уменьшаются в 
3-5 раз. В условиях полубесконечной диффузии по- 
лученные результаты совпадают с расчетами [1-3]. 

Для средних значений Н расчеты ВАК совпада- 
ют с известной теорией Де Вриза и Ван Далена [1] 
и работами других авторов [3,4]. Так же как и в ука- 
занных работах, нами была найдена область значе- 
ний Н, где параметры анодного пика (максимум 
пика и ширина пика на полувысоте) слабо зависят 
от величины Н. 

Замечено, что при Ж0,1 зависимость параме- 
тров анодного пика от Н более сложная, чем полу- 
ченная Де Вризом и Ван Даленом [1]. Получаемые 
пики в данных условиях более высокие и узкие, чем 
в [1]. Изучение области очень малых значений Н 
особенно актуально для развития метода инверси- 
онной вольтамперометрии в связи с общей тенден- 
цией к миниатюризации оборудования для прове- 
дения анализа. Полученные результаты в пределе 
согласуются с результатами теоретических и экспе- 
риментальных работ для микроэлектродов. 

Выводы 

1. Представлена физико-математическая модель 
обратимого электродного процесса для простой 
электрохимической реакции, протекающей при 
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линейном изменении потенциала на электродах 
ограниченного объема. 

2. Впервые предложены параметры Яи Ѳ 0 в каче- 
стве критериев подобия для плоских тонкоплё- 
ночных электродов, характеризующие массооб- 
мен на границе раздела фаз между электродом и 
раствором электролита. 

3. Теоретически показано, что величина Н учиты- 
вает соотношение скорости электрохимических 
реакций (прямой и обратной), определяемой 
скоростью изменения потенциала, и диффузии, 
обусловленной взаимодействием между моле- 
кулами в растворах (металлических, жидких) и 
квадратом толщины плёнки электрода (анод- 
ный процесс) или электролита (катодный про- 
цесс). Величина 0 О учитывает соотношение ско- 
ростей электрохимических реакций в зависимо- 
сти от равновесного потенциала. 

4. Показано, что при замене интеграла в точном 
уравнении вольтамперных кривых на производ- 
ную от концентрации по потенциалу время рас- 
чётов можно уменьшить на 2-3 порядка без су- 
щественной потери точности расчётов в обла- 
сти значений параметра Н от нуля до единицы, 
в практической области метода инверсионной 
вольтамперометрии . 

5. Выразив коэффициент диффузии в формуле 
для параметра Н через коэффициент вязкости, 
полученные соотношения можно применять не 
только для амальгамных, но и для твёрдых мо- 
дифицированных электродов. 
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